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L  ATTRACT* 

A  constant  discharge  coefficient  (C,j),  relating  flow  rate  to  the  square  root 
of  the  differential  pressure  across  a  primary  element,  is  a  desirable  characteristic 
of  any  flow  sensor  to  be  used  in  shipboard  automatic  control  systems.  Although  the 
simple,  square  edged  orifice  meets  this  criterion  for  turbulent  flows,  c£  does  not 
remain  constant  for  applications  in  the  laminar  and  transitional  regions  of  pipe 
Reynolds  number  (Rep).  In  a  quest  for  a  device  which  exhibits  both  simplicity  and 
a  constant  0^  at  low  Rep,  twelve  conic  entrance  orifice  plates  were  tested.  > 

The  range  of  orifice  to  pipe  diameter  ratios  (B)  was  0.1  to  0,5,  and  Rep 
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below  V^t  critical  Reynolds  number  for  pipe  flow  (~  2000).  Theoretical  considerations! 
showed  that  this  hump  is  attributable  to  fluid  viscosity  and  that  its  height  might  be 
lessened  by  geometrical  modifications.  Accordingly,  qualitative  arguments  are  pre¬ 
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R*p. 
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ABSTRACT 

A  constant  discharge  coefficient  (C^),  relating  flow  rate  tc 

the  square  root  of  the  differential  oressure  across  a  primary  element, 

is  a  desirable  characteristic  of  any  flow  sensor  to  be  used  in  shipboard 

automatic  control  systems.  Although  the  simple,  square  edged  orifice 

meets  this  criterion  for  turbulent  flows,  C.  does  not  remain  constant 

a 

for  applications  in  the  laminar  and  transitional  regions  of  pipe 
Reynolds  number  (ReD).  In  a  quest  for  a  device  which  exhibits  both 
simplicity  and  a  constant  at  low  Re^,  twelve  conic  entrance  orifice 
plates  were  tested. 

The  range  of  orifice  to  pipe  diameter  ratios  (jS)  was  0.1  to  0.5,  =ind 
ReD  ranged  from  40  to  50,000.  was  constant  within  +  1%  over  the 
^  specified  range  of  ReD  for  ft  -  0.1.  Bor  R  >  0.2,  Cd  was  constant  within 

+  2/6,  indicating  a  worsening  of  performance  with  increasing  ft. 

Graphs  of  versus  Re^  exhibited  a  ''hump"  which  consistently  occurred 
Just  below  the  critical  Reynolds  number  for  pipe  flow  (~  2000) .  Theoretical 
considerations  showed  that  this  hump  is  attributable  to  fluid  viscosity 
and  that  its  height  might  be  lessened  by  geometrical  modifications. 
Accordingly,  qualitative  arguments  are  presented  which  contend  that  a 
protruding  conic  edge  might  improve  performance  at  low  ReD. 

The  conic  e.  trance  orifice  is  to  be  preferred  to  the  square  edged 
orifice  where  low  values  of  ReD  are  encountered.  The  device  is  not 
unique,  however,  since  quadrant  edge  orifices  are  also  well  suited  for 
such  measurements. 

*  .v 
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SUMMARY  PAGE 

The  Problem 

A  constant  discharge  coefficient  (Cd),  relating  flow  rate  to  the 
square  root  of  the'  differential  pressure  across  a  primary  element, 
is  a  desirable  characteristic  of  any  flow  sensor  to  he  used  in 
shipboard  automatic  control  systems.  Although  the  simple,  square  edged 
orifice  meets  this  criterion  for  turbulent  flows,  Cd  does  not  remain 
constant  for  applications  in  the  laminar  and  transitional  regions  of 
pipe  Reynolds  number  (Re^) .  A  modified  orifice  for  which  Cd  remains 
constant  at  low  Rep  (i.e.,  into  the  laminar  region)  would  thereby 
contribute  to  Naval  flow  measuring  capabilities.  Accordingly,  it  is 
the  aim  of  this  project  to  investigate  the  performance  of  the  coziic 
entrance  orifice  as  a  flow  measuring  element  in  the  low  Rep  domain 
covering  laminar,  transitional  and  turbulent  pipe  flows. 

A  secondary  aim  is  to  decipher  which  modifications  to  this  type 
of  orifice  are  likely  to  result  in  improved  performance  over  the 
specified  Reynolds  number  range. 
findings 

While  the  conic  entrance  orifice  is  an  improvement  over  the  square 
edged  orifice  for  flow  measurement  in  the  low  Reynolds  number  domain,  it 
is  only  for  the  ana lie st  orifice-to-pipe  diameter  ratio  tested,  8  =  0.1, 
that  Cd  may  safely  be  assumed  constant  within  +  1%  maximum  variation. 

R>r  8  >  0.2,  the  discharge  coefficient  may  be  considered  constant  within 
+  2/6.  Further,  a  plot  cf  discharge  coefficient  versus  fi  with  less  than 
1%  standard  deviation  does  not  appear  feasible  for  the  specifial range  of 
Reynolds  number. 

Theoretical  considerations  suggest  that  modification  of  the  upstream 
surface  of  the  conic  entrance  orifice  plate,  producing  a  protruding  conic 
edge,  could  result  in  improved  performance  at  low  Reynolds  numbers. 

ii 
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teco™*>ndations 


The  conic  entrance  orifice  is  preferable  to  a  square  edged  orifice 

% 

wherever  pipe  Reynolds  numbers  below  2000  are  common.  The  device  is 

not  unique,  however,  since  the  quadrant  edge  orifice  shows  approximately 

» 

the  same  percentage  maximum  variation  in  discharge  coefficient  over  the 
specified  range  of  Reynolds  number. 

Further  investigations  should  involve  testing  of  modified  conic 
entrance  orifice  plates,  such  as  the  protruding  conic  edge  orifice 
proposed  in  this  report. 


iii 


i  V^«as>,a4^*WL,*  ^^'J*  it 


mavsecpfiiadiv  mojEcr  a-iooo 


NCHENCIATURE 

A  -  area 

Cd  -  discharge  coefficient 
D  -  pipe  diameter 
d  -  orifice  diameter 
b  -  3*88  flow  rate 
P  -  pressure 

Rep  -  pipe  Reynolis  number 

Re,  -  throat  Reynolds  number 
a 

v  -  fluid  Telocity 

flag* 

5  -  the  ratio  d/D 
p  -  fluid  mass  density 
Subscripts 

1  -  denotes  upstream  station 

2  -  denotes  downstream  station 
a  -  actual 

t  -  theoretical 
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INTRODUCTION 

\ 

The  comnion  square  edged  orifice,  consisting  of  a  concentric  hole  in 
a  circular  metal  plate,  is,  geometrically,  an  exceedingly  simple  device  for 
the  measurement  of  fluid  flow.  If  the  Reynolds  number  (Re^)  of  the  flow 
is  sufficiently  large,  then  application  of  Bernoulli's  equation  leads  to 
a  theoretical  expression  for  the  mass  flow  rate  which  is  equally  simplistic, 
and  is  directly  proportional  to  the  actual  mass  flow  rate.  The  "constant1' 
of  proportionality,  .  is  known  as  the  discharge  coefficient,  and  is 
actually  a  weak  function  of  the  Reynolds  number,  asymptotically  approaching 
a  constant  value  as  ReD  -  If  the  upstream  edge  of  an  ordinary  orifice 
is  beveled,  it  is  called  a  conic  entrance  oiifice.  Fbr  sufficiently  large 
Ref)  'rych  an  orifice  lends  itself  to  the  same  simplified  flow  measuring 
procedure  described  above  for  the  squared  edged  orifice.  However,  in  the 
low  Reynolds  number  domain,  v/hen  C<j  can  no  longer  safely  be  assumed 
independent  of  ReD,  the  dependence  of  on  ReD  is  different  for  the  square 
edged  and  coni^  entrance  orifices. 

Since  the  magnitude  of  the  deviation  of  from  its  mean  determines 
the  lower  lir-dt  of  Re^  for  which  the  Bernoulli  equation  approach  will 
lead  to  accurate  flow  measurement,  any  modification  to  an  ordinary  orifice 
which  results  in  a  smaller  variation  in  at  low  Reynolds  numbers  thereby 
results  in  an  orifice  whose  practical  range  of  applicability  is  increased 
(i.o.,  it  may  be  used  for  smaller  ReD).  The  primary  purpose  of  this 
investigation  is,  therefore,  to  experimentally  determine  the  values  of 
the  discharge  coefficient  for  various  conic  entrance  orifices  and  thereby 
assess  their  adequacy  for  flow  measurement  over  an  unusually  wide  range 
of  Reynolds  numbers,  from  the  laminar  region  into  the  turbulent  region. 
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A  secondary  objective  is  to  investigate,  analytically,  modifications 
to  the  conic  entrance  orifice  which  might  further  extend  the  range  of 
Re^  for  which  orifice  flow  measurement  is  practicable. 

The  conic  entrance  orifice  has  been  selected  bv  the  International 

t 

Standards  Organisation  (ISO),  Working  Group  IS0/TC30/WG1O,  as  a  recommended 
ISO  standard  when  flows  have  pipe  Reynolds  numbers  in  the  region  shown 
in  Table  I.  The  results  of  this  investigation,  in  addition  to  the  test 
results  of  similar  investigations  by  form?  leading  manufacturers  of  flow 
measurement  equipment,  will  be  made  available  to  the  Fluid  Meters  R- search 
Committee  of  the  American  Society  of  Mechanical  Engineers,  whir'u  wi1!  aid 
the  U.S.A.  committee  for  IS0AC30  to  establish  its  position  as  to  the  adequacy 
of  the  oonic  entrance  orifice  as  an  international  standard  for  low  Reynolds 
number  flow  measurements. 


TABLE  I  -  RANGE  OF  PIPE  REYNOLDS  NUMBERS  AS 
A  FUNCTION  OF  DIAMETER  RATIO  R 


ft 

Minimum 

Ren 

— - — - 

Maximum 

0.1 

40 

0.2 

40 

0.3 

60 

0.4 

120 

0.5 

260 

2 


1  t 

v.  * 


NAVSECPHIIADIV  H>OJSCT  A-1000 


SCME  THEORETICAL  BACKGROUND 

The  theoretical  equations  for  liquid  flow  through  a  conic 
orifice  are  normally  obtained  in  the  same  fashion  as  for  any  other 
differential  pressure  meter.  No  allowance  is  made  for  viscous  effects, 
a  flat  velocity  profile  is  assumed,  and  the  liquid  is  considered  in¬ 
compressible.  Using  this  model,  the  fundamental  equations  are  those 
of  mass  conservation, 

P\k  *  Pva4 

and  energy  conservation  at  constant  elevation: 

Pi-lpVi2  r  ?2+  2pVa2 

Denoting  the  ratio  of  the  orifice  diameter,  d,  to  the  pipe  diameter, 

D,  by  j5,  it  is  a  matter  of  simple  algebra  to  solve  the  above  equations 
for  the  theoretical  mass  flow  ra+e : 


N  *  pv^  - 


si 


If  the  actual  mass  flow  rate  is  denoted  by  1%,  then  the  discharge  coefficient 
is  defined  by: 

“a 


cd  *  nit 


Hence,  if  the  values  of  d,  D,  and  p  are  known,  measurement  of  the  difference 
between  upstream  pressure  (Px )  and  the  pressure  of  the  fluid  issuing  from 
the  orifice  (P2)  will  permit  determination  of  m^.  The  actual  flow  rate  is 
then  predictable  if  C^  is  known. 

Although  a  more  complete,  general  theory  of  fluid  flow  is  manifest 
in  the  Navier-Stokes  equations,  these  equations  have  never  been  solved 
exactly  for  flow  through  an  orifice.  The  preceeding,  simplified  method 
has  proven  most  practical  for  measuring  turbulent  flows,  for  is  then 
nearly  constant.  It  is  in  the  laminar,  or  low  Reynolds  number  domain  that 
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the  dependence  of  on  Re^  becomes  pronounced.  In  this  circumstance, 

one  must  know  the  value  of  ReD  in  order  to  predict  and,  hence,  ir^. 

But  m  must  be  known  in  order  to  calculate  Ren.  The  lamentable  situation 
a  u 

thus  arises  whereby  one  must  know  mfi  in  order  to  predict  m&.  With  this 
realization,  the  desirability  of  an  orifice  flowmeter  whose  discharge 
coefficient  e:chibits  only  slight  change  with  changes  in  Re^  in  the  laminar 
region  is  immediately  evident. 

EXPERIMENTAL  PROGRAM 
Test  Equipment 

The  geometry  of  a  conic  entrance  orifice  plate  is  depicted  by 
Figure  1.  The  labels  on  the  figure  identify  parameters  whose  specifications 
are  listed  in  Appendix  A.  Figure  2  is  a  photograph  of  the  upstream  face 
Of  one  of  the  conic  entrance  orifice  plates  used  in  the  experiments.  The 
conic  entrance  is  accented  by  the  light  reflected  from  it. 

Twelve  orifice  plates  were  tested  in  all.  Rive  were  manufactured 
by  Daniel  Industries,  Inc.,  and  seven  by  Taylor  Instrument  Division,  Sybron 
Corporation.  Two  inch  piping  (70"  upstream,  10"  downstream)  was  supplied  by 
each  manufacturer  for  use  with  -their  respective  orifices.  The  orifice-to-pipe 
diameter  ratio  ($)  had  nominal  values  of  0.1,  0.2,  0.3,  0.4,  and  0.5.  The 
dimensions  of  these  plates,  as  well  as  the  required  specifications,  are 
listed  in  Appendix  A.  It  should  be  noticed  that  the  mamfacturers  evidently 
found  certain  specifications  difficult  to  meet. 

Two  centering  dowels  proirudod  from  the  flar.gee  on  the  Daniel  piping 
so  that  the  outside  diameter  of  an  orifice  plate  (as  depicted  in  Figure  2), 
when  held  flush  with  the  dowels,  would  be  concentric  with  the  piping.  The 
Taylor  flanges  and  orifices  had  matching  pairs  of  holes  in  which  centering 
pins  (see  Figure  3)  were  to  be  placed  to  assure  concentricity.  A  paper 
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TGURE  3  -  CROSS  SECTION  DIAGRAM  OF  TAYLOR  FLANGE  ASSEMBLY 
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"by  Miller  and  Kneisel  joints  out  that  eccentricity  of  oil  tit  ■'  -dth 
piping  aay  effect  considerably  the  discharge  coefficient,  he-  :ns  necessitating 
these  precautions, 

in  order  to  attain  tty*  low  Reynolds  numbers  requisite  to  tnis 
investigation,  while  maintaining  differential  pressure  large  enough 
to  permit  accurate  determination  of  C^,  it  ras  necessary  So  employ  a 
relatively  viscous  fluid.  On  the  other  hand,  to  attain  the  highest 
Reynolds  numbers  called  for  (see  Table  I),  virile  avoiding  pressures  too 
high  to  measure  on  a  one-hundred  inch  mercury  manometer,  considerably  lower 
viscosities  were  necessary.  Preliminary  calculations  showed  that  water 
was  unsuitable  for  the  specified  range  of  ReQ  under  the  restrictions  on 
the  differential  pressure  range  (10"  H2<><  AP  <  100"  Hg).  Additional 
calculations  showed  that  the  specified  ranges  of  Re^  (subject  to  the 
restrictions  on  AP)  could  be  covered  through  the  use  of  two  oils  of 
different  viscosities:  Navy  Special  Fuel  Oil  for  low  Re^,  and  the  Navy's 
distillate  fuel  for  the  higher  values  of  Re^. 

Since  both  the  water  and  the  mercury  manometers  used  in  thes^  »sts 
required  ~  a  ter  billed  input  leads,  it  was  necessary  to  have  an  inter' 
between  the  metered  liquid  (oil)  and  water.  This  was  accomplished  by 
manufacturing  a  pair  of  oil-over-water  reservoirs  (see  Figure  4),  the  tops 
of  which  were  fed  from  upstream  and  downstream  corner  pressure  taps 
respectively  (see  Figure  3).  The  cross  section  areas  of  these  reservoirs 
were  made  large  enough  so  that  changes  in  manometer  water  levels  would  not 
appreciably  alter  the  height  of  the  interface  during  testing. 

Test  Procedure 

In  order  to  calculate  the  discharge  coefficient  accurate  iy.  it  was 
■^Superscripted  numbers  refer  to  references  listed  in  the  bibliography. 
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FIGURE  4  “  OIL/WATER  INTERFACE  RESERVOIRS 
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neoessaxy  to  obtain  precise  measurements  of  differential  pressure,  fluid 
temperature,  end  mass  flow  rate.  The  latter  parameter  was  Inferred 
indirectly  from  measurements  of  the  total  weight  flow,  the  total  time 
of  flow,  and  the  oil  specific  gravity  at  the  te^jerature  of  the  test.* 

The  differential  pressure  across  the  orifice  was  measured  with  one-hundred 
Inch,  u-tube  msnew*. .  r,  air  over  water  far  the  low  flow  rates,  and  water 
over  jarcury  for  tv  higher  flow  rates.  Fluid  temperature  was  measured 
using  iron-oonstantan  thermocouples  attached  to  the  surface  of  the  piping 
immediately  proceeding  the  orifice  and  about  65  Inches  upstream  (see 
Figure  5).  lagging  was  affixed  to  the  upstream  piping  (covering  both 
thermocouple  attachments)  so  that  transient  temperature  conditions  could 
be  detected  ty  the  accompanying  temperature  disparity  (i.e.,  at  low 
flow  rates  the  heated  oil  produced  a  temperature  rise  on  the  upstream 
thermocouple  considerably  in  advance  of  &he  temperature  rise  detected  near 
the  orifice).  Accurate  temperature  measurement  was  necessary  to  determine 
the  oil  viscosity,  and  hence  the  Reynolds  number,  accurately. 

The  accuracies  sought  in  the  measurement  of  primary  physical 
parameters  were  chosen  so  that  Cd  might  be  determined  within  +  1/4%, 
However,  due  to  fluctuations  in  the  flow;  the  necessity  to  measure  very 
email  differential  pressures;  large  response  times  of  temperature  and 
differential  pressure  elements  at  low  flow  rates;  and  unexpected  timing 
inaccuracies  for  very  fart  runs,  a  safer  estimate  for  the  mean  uncertainty 
in  the  determined  values  of  Cd  is  +  0.4%  (barring  nyst emetic  errors). 
Repeatability  checks  on  several  points  substantiated  this  estimate.  The 
uncertainty  in  ReQ  was  within  ±  10%. 

*  Oil  specific  gravity  was  measured  at  «p¥  and  calculated  forother 
temperatures  using  the  API  tables. 
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Test  Results 


All  twelve  orifice  plates  were  tested,  and  the  resultant  data 
were  used  to  calculate  discharge  coefficients,  C^,  and  pipe  Reynolds 
numbers,  Rep.  Dimensional  considerations  show  that  Cd  is  a  function 
only  of  Rep  and  the  plate  geometry.  The  results  of  this  investigation 
are  therefore  aptly  conveyed  by  a  graph  of  Cd  versus  Rep  for  each  orifice 
plate.  Such  graphs  are  given  in  Figures  6  -  17. 

Although  the  value  of  varies  with  it  is  characteristic  of 
all  conic  entrance  orifices  tested  that  Cd  has  an  absolute  maximum  some¬ 
where  in  the  laminar  region  (Rep  <  2000).  Although  this  investigation 
was  not  concerned  with  flow  at  very  high  Reynolds  numbers,  It  is  well 
■known  that  Cd  remains  nearly  constant  for  turbulent  flows.  It  appears 
as  though  the  ranges  of  Rep  listed  in  Table  I  are  such  that,  at  the  upper 
limit  of  Rep,  approaches  the  value  it  would  have  for  fully  developed 

turbulent  flow.  The  lower  limits  on  Re_  seem  to  have  been  chosen  so 

D 

as  to  exclude  values  of  C  that  are  lower  than  those  found  at  the  upper 

d 

limit.  The  height  of  the  "hump"  above  the  turbulent  value  of  Cj  may 
therefore  be  considered  as  a  measure  of  the  variation  in  Cd  over  the 
prescribed  range  of  Re^.  The  percentage  variation  in  Cd  is  then  gotten 
ty  dividing  this  variation  by  the  mean  value  of  in  the  turbulent  region. 
This  information  is  tabulated  below  (for  the  prescribed  range  of  Rep) . 


TABLE  I x  -  HRCL-.3TAGE  VARIATION  IN  Cd 

OVER  PRESCRIBED  RANGE  OF  Re„ 
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TABIZ  II  -  PERCENTAGE  VARIATION  IN  C a  OVER 


It  is  inters  sting  to  notice  that  Taylor  Plate  #7,  whose  discharge 


coefficient  exhibits  an  unusually  large  %  variation,  is  also  characterized 

by  a  throat  thickness  (dimension  "e"  in  Appendix  A.)  which  is  more  than 

fbur  times  the  value  called  for  try  specifications.  More  generally,  the 

pla+es  with  smaller  diameter  ratios  appear  to  show  less  deviation  in  C^. 

Fbr  practical  purposes,  when  a  conic  orifice  is  used  over  the  specified 

range  of  Re^,  should  he  considered  to  have  the  value  lying  midway 

between  the  extreme  values  of  C.,  thereby  having  an  uncertaintv  of  one- 

a 

half  of  the  $  variation  listed  above  (i.e.,  for  =  0.1,  Cj  might  be 
considered  accurate  within  ±1%).  It  should  also  he  noted  that,  for  the 

low  Reynolds  number  range  of  data  compiled  'ey  the  American  Society  of 

2 

Mechanical  Engineers  on  square  edged  orifice  plates,  there  is  a  greater 
%  variation  in  than  was  found  for  the  conic  entrance  orifices  described 
herein  (for  the  same  Reynolds  number  range). 
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Although  all  of  the  teat  results  are  embodied  in  Figures  6-17, 
it  is  nonetheless  informative  to  plot  the  dir charge  coefficient  versus 

t 

the  Reynolds  number  based  on  the  orifice  diameter,  Re^  (throat  Reynolds 
number).  Figure  18  is  such  a  plot,  ana  the  similarity  of  the  curves  for 
the  various  orifices  is  quite  apparent.  In  particular,  the  value  of 
Re^  for  which  is  a  jsaxirasr  is  essentially  the  same  for  all  the  orifices 
tested,  regardless  of  the  value  of  ft.  The  implication  is  that  Red  is 
a  physically  more  meaningful  dimensionless  group  against  which  to  plot 
Cd.  This  follows  if  one  interprets  the  change  in  sign  of  the  slope 
of  the  graph  as  signifying  the  transition  from  laminar  to  turbulent 
flow  in  the  orifice.  With  this  assumption.  Figure  18  indicates  that 
Red,  rather  than  ReD,  more  completely  cliaracterizes  thp  nature  of  the  flow 
through  the  orifice,  with  the  effect  of  R  being  merely  to  lower  or  raise 
the  characteristic  curve. 


A  more  exact  treatment  of  energy  conservation  than  the  usual 
application  of  Bernoulli's  equation  may  be  used  to  estimate  the 
dependence  of  the  discharge  coefficient  on  ft.  Appendix  C  should  be 
consulted  for  such  an  analysis. 
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ANALYSIS 

Although  simplicity  -characterizes  the  morphic  description  of  an 
orifice  plate,  its  dynamic  behavior  is  nothing  less  than  complicated. 

In  fact,  an  exact  solution  for  the  flow  of  a  real  fluid  through  an 

orifice  has  never  been  attained.  TO  accomplish  this,  the  Navier-Stokes 

r 

equations  would  have  to  be  solved  for  the  appropriate  boundary  conditions 
(i.e.,  the  velocity  vanishes  on  all  solid  surfaces).  However,  the 
nonlinearity  of  these  equations,  as  well  as  the  discontinuous  boundary 
conditions  has  made  their  solution  unfeasible .  A  more  general  and  more 
qualitative  approach  must  be  employee,  if  the  performance  of  the  conic 
orifice  Is  to  be  ''understood".  The  most  reasonable  aim  of  such  an  under¬ 
taking  is  to  decipher  precisely  which  modifications  to  the  plate  geometry 
will  result  in  a  lesser  variation  of  in  the  low  Reynolds  number  domain. 

Dimensional  analysis  of  the  parameters  which  characterize  incompressible 
flow  through  a  conic  orifice  shows  that  a  complete  description  of  any  test 
with  a  conic  orifice  is  expressible  in  terms  of  the  dimensionless  numbers: 

C^,  Rep,  /?,  F,  ~  and  ^  (the  last  three  symbols  refer  to  Figure  1).  Further, 
if  the  viscosity  is  neglected,  then  Re_  is  no  longer  necessary.  The  important 
conclusion  is  that,  for  an  ideal  (inviscid)  fluid,  is  a  function  of  the 
plate  geometry  and  P  alone.  Since  Cd  would  be  constant  for  all  flow  rates 
were  the  fluid  inviscid,  it  must  be  the  viscosity  which  causes  the 
characteristic  "bumps"  in  Figures  6  -  17. 

It  is  instructive  to  solve  the  problem  of  ideal  fluid  flow  (potential 
flow)  through  an  ordinary,  square  edged  orifice.  This  can  be  done  for  the 
case  of  vanishing  /?,  with  the  added  assumption  of  uniform  mass  flux 
through  the  plane  of  the  orifice.  (These  assumptions  greatly  facilitate 
solution.)  The  problem  is  worked  out  in  Appendix  B.  The  solution  thus 
obtained  predicts  a  discharge  coefficient  of  approximately  0.65.  Under 
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the  assumptions  made,  the  direction  of  the  fluid  velocity  at  the 

orifice  edge  necessarily  makes  an  angle  of  about  45°  with  the  pipe 

axis.  To  assess  the  validity  of  this  model,  the  experimental  results 

must  be  consulted.  The  average  value  of  the  discharge  coefficient 

for  slightly  turbulent  flow  (RCp  ~  2000)  through  an  orifice  with 

p  =  0.1  is  about  0.612.  Although  the  prediction  given  above  is 

about  7%  too  high,  it  should  be  noted  that  the  analysis  in  Appendix  B 

shows  that  the  exact  prediction  is  actually  less  than  0.65  (perhaps 

o 

by  several  per  cent).  It  is  also-  interesting  to  note  that  45  is  the 
optimum  edge  bevel  angle  for  conic  entrance  orifices  with  small  P. 

Based  on  the  above  results,  it  is  clear  that  the  analysis  and  assumptions 
made  form  a  reasonably  accurate  description  of  the  actual  flow  through 
a  square  edged  orifice  wi&h  small  diameter  ratio  («)  in  the  region  of 
low  turbulence  (Rep  ~  200G). 

Although  tl '  above  results  are  gratifying,  they  do  not  point 

directly  to  e  means  of  improving  the  conic  entrance  orifice.  Still, 

it  should  be  kept-  in  mind  that  the  more  closely  the  test  conditions 

app..  ,ach  the  ideal,  the  more  nearly  might  be  expected  to  remain 

constant.  Before  discussing  plate  modifications  which  might  be  expected 

to  result  in  mere  nearly  ideal  boundary  conditions,  the  reason  for  the 

improved  performance  of  ‘..he  conic  entrance  orifice  as  compared  to  f 

square  edged  orifice  should  be  elaborated. 

The  exact  solution  of  the  Navier-Stokes  equations  for  the  steady 

flow  of  a  real,  incompressible  fluid  (liquid)  along  a  converging,  plane- 

3 

walled  channel  was  obtained  by  G.  Hamel  in  1916  .  Although  the  solution  is 
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very  complicated  mathematically,  qualitatively  the  results  are  that,  for 
ReE  much  larger  than  1  (perliaps  several  hundred),  the  flow  is  essentially 
lilce  that  of  an  ideal  fluid  (i  .e . ,  potential  flow) .  Although  this  solution 
is  for  the  two  dimensional  case,  it  is  very  reasonable  to  assume  that  its 
essential  characteristics  would  be  found  in  the  exact  solution  for  convev-'':nrr 
flow  in  a  cone.  Thus,  it  would  be  expected  that  the  introduction  of  n  c^- 
entrance  to  a  square  edged  orifice  could  be  considered  as  a  perturbation  ■-■■■ 
the  flow  whose  effect  is  to  bring  the  flow  pattern  into  closer  correspond^  ee 
to  the  potential  flow  solution.  It  should  be  remembered  that  this  reason^  ~ 
is  valid  only  for  ReD  »  1 . 

Mathematically  speaking,  the  only  requirement  for  potential  flow 
is  that  the  curl  cf  the  velocity  vanish  throughout  the  region  under 
consideration.  Physically,  non-zero  curl  is  created  in  regions  of  fluid 
shear  stresses  which  are  attributable  to  non-zero  viscosity  plus  the  no-slip 
boundary  condition  (i.e,  v  =  0  on  all  surfaces).  Since  the  fluid  viscosity 
Is  something  that  cannot  be  eliminated,  it  is  only  through  the  boundary 
conditions  that  the  curl  of  the  velocity  might  be  minimized.  R>r  practical 
values  of  Rep,  zero  curl  is  never  a  good  approximation  downstream  from  an 
orifice  plate.  Therefore,  only  modifications  to  the  upstream  boundary 
conditions  will  be  considered  here. 

Examination  of  Figures  6-17  shews  that  increasing  p  has  an 
undesirable  effect  on  the  behavior  of  C^.  Since  the  magnitude  of  p  is 
indicative  of  the  relative  proximity  of  the  pipe  wall  to  the  orifice,  it 
it  reasonable  to  surmise  that  the  nearness  of  bounding  surfaces  to  the 
orifice  entrance  enhances  non-ideal  performance .  It  appears  that  a  6tep 
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•toward  ideal  performance  would  involve  a  decrease  in  the  surface  area 
near  the  conic  entrance.  Aside  from  the  pipe  walls,  the -only  other 
surface  near  the  conic  entrance  is  the  upstream  face  of  the  orifice 
plate  itself.  This  may  therefore  he  accomplished  in  a  straightforward 

i 

manner  iy  simply  recessing  that  portion  of  the  face  of  a  conic  entrance 
orifice  (Figure  1)  which  is  nearest  the  orifice  entrance.  A  conic 
entrance  orifice  plate  thusly  modified  might  he  called  a  ''protruding 
conic  edge"  orifice  plate,  a  diagram  of  which  is  Figure  19.  It  should 
he  emphasized  that,  to  the  best  of  the  author's  knowledge,  such  an 
orifice  does  not  presently  exist,  and  Figure  19  is  merely  one  possible 
design  which  seems  worthy  of  investigation.  The  central  ideal  behind 
this  type  of  modification  is  to  decrease  the  influence  that  the  viscous 
effects  associated  with  the  plate  surface  have  on  the  orifice  inlet 
velocity  profile . 

It  should  also  he  mentioned  that,  although  theoretical  considerations 
might  lead  to  an  extended  range  of  applicability  of  oi*ifice  flow  measurement, 
they  also  lead  to  the  conclusion  that  orifice  flow  measurement  cannot  be 
practicable  in  the  limit  of  very  small  Reynolds  numbers  (Re^  <  1).  Algebraic 
manipulation  of  the  exact  solution  for  slow,  viscous  flew  through  a  slit** 
shows  that  must  vanish  in  the  limit  of  very  slow  flow  (v  0).  In  other 
words,  although  it  is  reasonable  to  attempt  to  eliminate  the  characteristic 
hump  in  vs.  Re 3,  it  is  impossible  to  eliminate  the  mono tonic  decline 
in  for  values  of  Re^  to  the  left  of  the  hump. 
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CONCLUSIONS 


While  the  conic  entrance  orifice  is  an  improvement ' over  the  square 

edged  orifice  for  flow  measurement  in  the  low  Reynolds  number  domain 

(Table  I),  it  is  only  for  the  smallest  diameter  ratio  tested,  8  ~  0.1, 

that  may  1?e  expected  to  remain  constant  within  +  1%.  Further  comparison 

of  the  data  obtained  in  this  investigation  with  data  obtained  in  other 

investigations  indicates  that  a  plot  of  C^  (average)  vs.  R  with  less 

than  1%  standard  deviation  (as  sought  by  the  International  Standards 

Organization)  is  not  feasible  for  the  specified  Reynolds  number  range. 

Also,  the  conic  entrance  orifice,  is  not  unique,  since  the  quadrant  edge 

orifice  shows  approximately  the  same  percentage  maximum  variation  in 

* 

discharge  coefficient  over  the  specified  range  of  Reynolds  number0. 

For  Reynolds  numbers  above  800,  Cd  for  the  quadrant  edge  orifice  actually 
exhibits  a  lesser  percentage  variation. 

There  is  reason  to  believe  that  further  modification  of  the  conic 
entrance  orifice  plate,  such  as  the  protruding  edge  design  in  Figure  18, 
might  result  in  improved  performance.  Strict  specifications  would  he 
necessary  for  the  manufacturer  of  any  modified  plates.  In  particular, 
conformance  to  the  specification  for  edge  thickness  (parameter  "e") 
was  found  to  he  entirely  necessary  for  the  tested  orifices  (i.e.,  Taylor 
plate  #7,  R  =  0.2,  had  an  edge  thickness  of  about  four  times  the  specified 
value,  and  its  plot  of  vs.  Re^  was  anomalous  compared  to  the  other  two 
plates  having  8  -  0.2).  In  view  of  the  difficulties  manufacturers  had 
meeting  the  required  specifications  for  the  test  plates,  it  appears  as 
though  mairofacturing  problems  could  limit  practical  use  of  the  conic 
entrance  orifice. 

It  should  be  noted  that  conversion  from  the  Navy's  diesel  fuel  oil 
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to  Navy  distillate  will  lower  the  range  of  Reynolds  numbers  encountered 
in  certain  shipboard  systems.  This  development  could  give  added  value 
to  the  conic  entrance  orifice,  which  is  well  suited  to  low  Reynolds 
number  application. 

t 

RECOMMENDATIONS 

j..  The  conic  entrance  orifice  should  be  considered  on  a  par  with 
the  quadrant  edge  orifice  for  shipboard  flow  measurement  at  low 
Reynolds  numbers.  Practical  considerations,  such  as  cost  and  manu¬ 
facturers  ability  to  meet  specifications,  should  dictate  whether  the 
conic  entrance  orifice  should  he  used  instead  of  other  (low  Reynolds 
number)  orifices. 

2.  Where  flow  rates  are  encountered  which  considerably  overlap  the 
regions  listed  in  Table  I,  the  conic  entrance  orifice  is  to  be  preferred 
to  a  square  edged  orifice. 

3.  Further  investigations  aimed  at  low  Reynolds  number,  orifice  flow 
measurement  should  study  protruding  conic  edge  designs  such  as  that 
shown  in  Figure  19. 
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APPENDIX  A 

SPECIFIED  AND  MEASURED  VALUES  OF 
CONIC  ENTRANCE  ORIFICE  PIATE  DIMENSIONS 


I 


m 

m 


Ncvainal 

Beta 
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DIMENSIONAL  COMPARISONS  (INCHES) 


Taylor 

Plate 


.1186 

43P47‘ 

.0568 

.0152 

.6214 

1.2423 


.2997 


.1186 

39°28» 

.0800 

.0205 

.8271 

1.6534 

.3989 


Required 

Spec. 


J&L 

.126 

39?  +  1.5° 

.080 

.017 

.8275 

.8265 

IMk 

1.654 


Daniel 

Plate 


.1274 


.0535 

.0165 

.6196 

1.2426 

.3003 
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SIMPLIFIED  SGDTTIOfl  PCS.  POTE’.’I'IAL 
FLOW  THECUGi  AN  ORIFICE 
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NOMENCLATURE 

A  -  area 

a  -  radius  ox*  orifice 
Cjj  -  discharge  coefficient 
F  -  hypergeometric  function 
f  -  a  function  of  k 
G  -  a  function  of  n 
Jn  -  Bessel  function  of  order  n 
k  -  integration  variable 
P  -  pressure 
v  -  fluid  velocity 
2  -  coordinate  along  symetry  axis 
Greek 

-  ratio  of  orifice  diameter  to  pipe  diameter 
y  -  a  substitution  for  the  power  series  in  £ 

T  -  the  gamma  function 

Q  -  angle  of  converge  of  the  flow  at  the  orifice  edge 
£  -  substitution  .for  { 

0  -  radial  cylindrical  coordinate 
pm  -  mass  density 
<p  -  the  velocity  potential 
Subscripts 

o  -  at  the  plane  of  the  orifice 
1  -  at  the  plane  of  the  vena  contracts 
r  -  in  the  radial  direction 
t  -  total 

z  -  in  the  axial  direction 
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The  analysis  that  follows  treu.cs  x^e  flow  of  an  incompressible, 
lnvi3cii  fluid  through  a  common  orifice.  Only  the  upstream  velocity 
field  is  treated,  since  the  condition  of  zero  curl,  which  is  necessary 
for  potential  flow,  is  almost  never  a  reasonable  approximation  down¬ 
stream  of  an  orifice.  The  condition  of  incompressible  flow  (tx  very 
good  assumption  for?  liquids)  is  expressible  mathematically  as 

v  •  v  =  0 

Tho  condition  of  zero  curl  permits  the  expression  of  the  velocity 
as  the  gradient  of  a  velocity  potential: 

That  v  x  v  =  0 

permits  the  velocity  to  be  expressed  as 

v  =  Vo 

The  equation  which  the  potential  must  satisfy  throughout  the  upstream 
volume  is  therefore  Laplace’s  equation; 

V8<p  -  0 

Because  the  problem  has  cylindrical  syrame+zy,  a  general  mathematical 

B1 

expression  may  be  used  for  the  potential  in  the  semi-infinite  domain  : 

■n  (o,  z)  «  fdk  f(k)  e  k  "  J  (kp) 

o  Equation  (Bl) 

where  the  form  of  f(k)  depends  on  the  boundary  condi  liens .  If  the 

orifice  has  radius  a,  then  an  obvious  boundary  condition  is  that  the 

component  of  velocity  normal  to  the  plate  vanish  at  the  plate: 

vz  (o,0)  =  ~  0  for  a  <  p  <  ® 

To  completely  determine  the  solution  uhroughoui  the  region,  the  boundary 


B2 
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\VSfflyssawv  r -  .«- 


i  V.,'* 


condition©  must  be  specified  over  the  entire  bounding  surface.  Since 
the  geometry  doed  not  dictate  the  boundary  conditions  at  the  orifice . 
we  assume  that  the  efflux  of  fluid  is  uniform  over  this  surface: 

t 

vz  {p, 0)  =  =  v  for  0  <  o  ^  a 

*  dz  *o 

Noticing  that 

(o,0)  =J  dX  X  f(X)  Jc(Xo) 

the  boundary  conditions  which  determine  f(k)  are  seen  to  be: 

/* ’“'<*>  J0(to,)-{>  °<o<- 

To  find  f (k)  from  these  conditions  it  is  necessary  to  employ  the 
B2 

Hankel  transform  ,  which  states  that  if 


G(o)  =  ;  dk  k  f(k)  J0  (kp) 


then  the  inversion  formula  is: 


f(k)  =j  do  oG(n)  J0(ko) 


For  our  problem  we  have 


so  that 


0(r.)  -  (-%  0  *  0  *  " 
‘.0  a  <  n  <  °» 


f  (k)  “  -  [  do  P'vz  )J0(kp) 


which  may  be  evaluated  u  7ing  tne  property  of  Bessel  fmctions  that 


d  UnJnfx)'J  =  xnJn_1(x) 
dx 


B3 


\  i>3  " 
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! 


{  '> 
Vj>' 


-va 


v(p, o) « ^  +  f  e  *  ^  f8  +  igg  *lo+ . . .  i 

2ff  L  -  -  32  e  ?-152  h  73723  r  7372300  J 

Recalling  that  =  -,  and  ^  us  have : 

Sp  a'  ;  do  SA  5n 

vr(p,0)  =  [$  +  .375J3  +  ,2344^  +  .1709£7  +  .1346a3  Equation  (B2) 


Letting  y  represent  the  serf.es  terns  in  the  brackets,  the  velocity 
profile  on  the  plane  of  the  orifice  may  be  expressed  as: 


vo  “  vzq  (2  ~2-r) 


where  s’  and  'v  are  unit  vectors  in  the  axial  and  radial  directions 
respectively. 

Since  we  assume  potential  flow,  the  total  pressure  is  a  constant 
throughout  the  fluid,  and  depends  on  the  usual  (static)  pressure,  P, 
and  velocity  as  follows: 

pt  *  f  *  I  "m  * 

The  statement  that  this  is  a  constant  is  equivalent  to  stating 
that  energy  is  conserved.  The  remaining  conservation  lews  which  pertain 
bt°  momentum  con  Nervation  and  mass  conservation.  We  apply  these  laws 
between  the  plane  of  the  orifice  and  the  plane  of  the  vena  contracts, 
with  the  added  assumptions  that  there  is  no  flow  through  the  radial 
walls  of  our  control  volume  (a  cylinder)  and  that  tire  pressure  in  the 
P-len©  of  the  vena  contracts  is  a  constant  and  is  in  fact  the  pressure 
which  is  measured  in  a  real  orifice  test.  Assuror g  further  that  the 
velocity  pro.'.'ile  in  the  vena  ccntracta  is  flat,  the  area  of  the  vena 


n? 
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contracta  may  be  calculated .  The  ratio  of  this  area  to  the  orifice 
area  will  then  be  considered  as  our  predicted  value  of  the  discharge 
coefficient. 

Pbr  the  general  expression  of  momentum  conservation  page  13  of 
reference  B4  may  be  consulted.  For  this  problem,  we  denote  the  orifice 
plane  by  o  and  the  plane  of  the  vena  contracta  by  1.  Mamontun  con¬ 
servation  becomes: 

pa  1  pb  .a 

0J  [pt  +  "  s  2ffod"  V-[?t  f  277060  l JP'.2^ 

where  the  radius  of  the  vena  contracta  is  b.  Letting  A  denote  area 
we  get 

AqP^'zo  f.l 

'  ~J  f  =  -  (\  *  VV  *0  *  (?t  *  VI  >  4  *  h  («o  -  4  ) 


Before  going  further,  the  integral  must  be  evaluated.  Recalling  the 
definition  of  y  (y  =  £  *  .375  f  +  .2344 f6  +  .1709*7  +  .1346^),  the 
integration  is  cumbersome ,  but  straightforward.  The  result  is: 


rl 

y3fdt  =  0.5552 
oJ 

Substituting,  the  principle  ecuation  becomes 

-  Lif21  4,  =  +  Vv  Ao  *  <pt  *  ipm\  >  4  +  P.  (V  -  4  ) 

Substituting  P,  =  -  ir^v^  and  cancelling  terms  gives; 

-  0.1388  vz3  =  -  0.5  v*^  +  0.5  A  +  0.5  v*  (A,  -  A0) 


0.2776 


(Zs>?s  /Zatf  -/  a.  h  , ,  /jL\ 

W'  (-V 


W’f  v„  vvnv\v,*r„r- 
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recalllrig  that  continuity  reauires  v^A  =  v  Ax,  or  (Vgc>V)  =  (**■},  we 

°  “P-  'VZj/  \2^/ 

obtain,  on  rearranging  terms: 


/•A a\2  /At  \ 

(7-  -'2.768:5  (2.;+  1.3843  =  0 
*iL%/  >iu/ 


,  fh  \ 


V 

Since  the  ratio  (Aj  /Aq)  is  oust  the  theoretical  correction  to  the 
mass  flovr  rate  predicted  using  Bernoulli's  equation,  the  root  of  this 
equation  is  the  theoretical  discharge  coefficient: 


~  2,.768>_-  J(2.7685T-  4(1.3843) 


or 


Cd  =  *65 


This  result  should  be  compared  to  the  average  experimental  value 

of  Cj  with  R  -  0.1.  which  is  approximately  0.61.  It  should  be  noted 

that  retention  of  higher  order  terms  in  the  previously  obtained  (hyper- 

geometric)  series  would  result  in  a  smaller  predicted  value  for  C, 

d* 

Because  of  the  poor  convergence  of  this  series  near  o  =  a,  the  exact 

solution  (retaining  all  terms)  is  probably  a  few  percent  lower  than 

the  approximation  given  here.  It  is  therefore  quite  clear  that  the 

enaiyais  given  above  is  a  reasonably  accurate  description  of  the  actual 

flow  through  a  square  edged  orifice  with  small  diameter  ratio  (R)  in  the 

region  of  low  turbulence  (Red  ~  2000). 

The  angle  of  convergence  of  the  fluid  velocity  at  the  edge  of  the 

orifice  may  be  estimated  by  substituting  ?  =  1  into  equation  (B2) : 

v  c  v 
r  o 

so  that  the  angle  of  convergence  is  approximately 
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.  »  -  *»  X%)  * 

which  happens  to  he  the  optimum  bevel  angle  for  the  conic  entrance 
orifice  (with  Q  ~  0.1). r- 
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APPENDIX  C 

CORK}.,  HON  TO  DISCHARGE  COEFFICIENT 
DUE  TO  LAMINAR  UPSTREAM  VELOCITY  FROFIIE 
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In  the  analysis  tliat  follows  it  will  he  shown  that  the  principles 
of  energy  and  mass  conservation  can  be  used  to  derive  a  more  exact 
expression  for  the  flow  rate  through  an  orifice  than  that  noimally 
obtained  using  the  Bernoulli  equation  for  one- dimensional  flow.  Viscous 
energy  dissipation  will  be  neglected. 

If  e  denotes  the  internal  energy  per  unit  mass,  then  the  energy 
per  unit  volume  of  the  fluid  is : 

ipv3  +  pc 

The  time  rate  of  change  of  this  energy  density  is  simply 

dt  (ip^  +  pO 

Using  the  equation  of  mass  continuity, 

||  •  pv  =  ° 

along  with  Euler’s  equation  for  inviscid  flow, 

|f  ■=-(?.?)  v -i?p 

P 

we  find  that 

St  <*"**>  *  +  p5  •  at 

«  -4—4-4  —4  r*  — f 

-  -  £irV-  pv  -  v  •  vp  -  pv  •  (v  •  V)  v 

=  -  Vv-  pv  -  V  .  VP  -^V  . 

The  thermodynamic  definition  of  the  differential  change  in  enthalpy 
(w)  is: 

dw  =  Tds  +  (i^dP 

where  e  is  the  specific  entropy.  We  may  therefore  write 

Cl 


aagBaBgaaaaaaaa^a^aaa  - 


-.S'^'V'T^',  *%T.rt«5^« 
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^P  -  -  p$6 

which,  upon  substitution  in  the  expression  for  the  time  rate  of 
change  of  kinetic  energy,  yields 


at  (i  pv3  )  =  -lb3?  *  ov  -  pv  •  ^(iv2  +  w)  +  pTv  •  vs  Equation  (Cl) 
We  now  consider  the  term-^pf)  "by  writing  the  thermodynamic  relation 


at 

d€  *  Tds  -  Pd(-) 


Tds  + 


■pj 

(^)^P 


and  recalling  that  w  =  c  +  E ,  so  that 

P 


d(pe)  =  cdp  +  pdc 

=  (w  -  ^ )  dp  +  p(Tds  +  ^  dp) 

=  wdp  +  pTds 

which  means 

?lpfcl  =  *  &A  +  Dx  2L£ 
dt  dt  p  at 

In  the  absence  of  heat  conduction  and  viscous  energy  dissipation,  the 
entropy  per  unit  mass  of  a  fluid  "particle "  will  remain  constant  in  time: 

dt'lt  ♦*•*•>*0 

or 

=  -v  •  Vs 

at 

likewise,  the  continuity  equation  for  mass  may  he  written: 

It  •  Pr 

the  expression  for  djflf )  becomes,  with  the  above  substitutions. 


-£*r&  QEaaag 


mj^smaaMOassk 


C2 


r^nv^  »^^v^.«%»*^^(rS,4r^5^^ 


■Wf/WK'flp-W--  A-J^’  4- 
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Equation  (C2) 


d  fcg)  =  -*?  .  0V  -  pTv  •  Vs 

o*t 

The  time  rate  of  change  for  the  total  energy  per  unit  volume  may 
therefore  he  written  (using  equations  (Cl)  and  (C2) 

=  -  (*v2  +  w)$  .  ft^  -  ov  .  vdv2  +  w) 
or,  finally,  using  the  relation  V  •  f g  =  £V  •  g  +  g  •  Vf, 

a(jp^.tne)  =  i?.  [pV(iv=  +  »)] 

This  equation  nay  be  applied  to  steady  state  orifice  flow 

~  Oj  by  integrating  over  some  volume  and  transforming  this 
volume  integral  to  a  surface  integral  using  Gauss*  divergence  theorem: 

0  =  jv-  [ov(iv2+  w)]dV  =  fovt^v2  +  w)>  dA 
volume  enclosing 

surface 

let  At  he  a  plane  surface  normal  to  the  pipe  axis,  and  sufficiently 
for  upstream  of  the  orifice  that  v  is  parallel  to  the  pipe  axis,  and  let 
fe  he  the  surface  in  the  plane  of  the  orifice.  Using  the  fact  that 
v  =  0  on  the  pipe  surface  and  orifice  surface,  the  above  equation 
becomes 

[  (io^3  +  ovw)dAi  =  j  ovfiv2  +  *)  •  34 
pipe  cross-section  orifice 

Since  the  aim  of  this  analysis  is  to  correct  for  the  upstream 


velocity  profile,  the  same  simplifying  (but  in? c curate)  assumption  will 
be  made  about  the  flow  at  the  orifice  as  is  normally  made  in  the  application 
of  Bernoulli's  equation:  v  is  parallel  to  the  pipe  axis  at  the  orifice. 
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(Actually,  this  assumption  is  more  correctly  applied  to  the  vena 
* 

contracts . )  All  that  remains  is  to  assume  velocity  profiles  at  4 

and  Aq  ,  and  the  above  equations  may  he  integrated.  To  keep  the 

/ 

calculations  relatively  simple  yet  accurate,  consider  flow  for  which 
Re^  >  2000,  yet  Rep  <  2000.  In  this  region  laminar  flow  will  he 
encountered  at  4  ,  while  turbulence  will  exist  at  4 .  The  velocity 
profile  for  laminar  pipe  flow  is  well  known  to  he  a  paraboloid,  while 
a  flat  velocity  profile  has  been  found  to  be  a  reasonable  approximation 
for  turbulent  flow.  In  other  words 

at  v  =  2Vx(l  -§) 

at  4  :  v  =  v, 

where  V  means  average  velocity,  and  use  has  been  made  of  the  fact 
that  the  velocity  along  the  pipe  axis  is  twice  the  average  velocity 
for  laminar  pipe  flow.  The  central  equation  now  becomes : 

^40^(1  -  ({?/)  2rrrdr  +  pvxi*4  =  pv24  (iv|  +  ) 

0  J 

which  upon  integration  gives 

_  2  _  ,_2 

ovx4  (vx  +  Wi  )  *  dv24  (iv2  +  ^  ) 

mass  continuity  requires  pvx  4  =  pv24,  hence 
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u 


p 

the  enthalpy  is  w  =  **  +  €,  so  that 

k 

2  p  i  2  P 

x  +  +  €i  =  i?3  +  a.  +  €3 

P  p 

The  change  in  the  internal  energy  is 

/ 

dc  =  Tds  -  ^  dp 
Pa 

However,  the  fluid  is  being  considered  incompressible  and  isentropic. 
so  that 

ds  =  dp  =  dc  =  0 

which  means  that  s,  =  £a .  The  equation  of  energy  conservation 
becomes : 

2  2 
nv i  +  Px  =  iov2  +  P3 

Apply  the  continuity  equation,  nv^  ~  ovn  A, , 


vi  =  %  ~  =  V?(^')=  ^Vj» 


and  substituting  this  for  v1 , 


p7a(T  +  P,  =  20^2  +  P2 


solving  for  v2. 


_  rmzxr 

VS  =V  n(l  -  2/f) 

The  theoretical  mass  flow  rate  is  therefore 

fzotir  -  Ta ) 

m*  =  nVaAe  =  7T(f  V  (i  I  ) 

Recalling  the  theoretical  expression  obtained  on  page  3  of  the  text, 
and  defining  a  new  discharge  coefficient,  C^,  we  have 


/ 

>1000 


C5 
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t 


Or 


%  *  c<pt  =  CK 


£d  =  mt  _  /  i  -  2(F 

c,  7  V  i  -  V 

C  6 


This  correction  factor  is  tabulated  below: 

_flL 


0.1 

0.2 

0.3 

0.4 

0.5 


0.9999 

0.9992 

0.9959 

0.9868 

0.9661 


The  validity  of  the  preceding  analysis  is  brought  to  light  by 
a  table  comparing  the  newly  defined  discharge  coefficient  to  the 
unual  discharge  coefficient  in  the  region  which  this  analysis  treats, 
Bed  -  2000.  (This  is  where  the  maximum  of  Cd  occurs,  and,  re  is 


hypothesized,  where  the  transition  to  turbulent  orifice  flow  begins.) 

» 


JL 

Cdfmaxl 

Measured  bevel 

0.1 

.745 

.745 

45°05* 

0.2 

.752 

.751 

45g08» 

0.3 

.751 

.748 

44q24’ 

0.4 

.782 

.772 

40  294 

0.5 

.815 

.787 

32°10» 

These  are  nominal  p  values,  and  average  values  of  and  F 

for  the  plates  tested.  For  the  first  three  plates,  the  bevel  angle  is 
approximately  the  same,  and  Cg  is  the  same  (within  better  than  1%) .  For 


the  last  two  p  values  the  decreasing  angle  of  convergence  apparently  results 

i 

in  a  larger  value  for  C^. 

In  summary,  it  appears  as  though  the  increased  height  of  the  hump 
in  Cfl  for  larger  p  values  is  due  to  the  fact  that  the  usual  Bernoulli 


C6 
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equation  approach  does  not  tabs  account  of  the  true  upstream  velocity- 
profile.  It  is  expected  that  if  F  were  the  same  for  all  values  of  $, 

t 

the  curves  of  versus  Re^  wouSd  he  nearly  coincident  for  every  plate, 
provided  the  correct  velocity  profile  were  used. 

r 
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